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ABSTRACT 

We investigate star formation rate (SFR) calibrations in light of recent developments in the modeling of 
stellar rotation. Using new non-rotating and rotating stellar tracks of Ekstrom et al. (2012), we study the 
integrated properties of synthetic stellar populations and find that the UV to SFR calibration for the rotating 
stellar population is 30% smaller than for the non-rotating stellar population, and 40% smaller for the Ha 
to SFR calibration. These reductions translate to smaller SFR estimates made from observed UV and Ha 
luminosities. Using the UV and Ha fluxes of a sample of ~ 300 local galaxies, we derive a total (i.e., sky- 
coverage corrected) SFR within 11 Mpc of 120-170 Mq yr^^ and 80-130 Moyr^^ for the non-rotating and 
rotating estimators, respectively. Independently, the number of core-collapse supemovae discovered in the same 
volume requires a total SFR of 270lgJ" M© yr~^, suggesting a mild tension with the SFR estimates made with 
rotating calibrations. More generally, when compared with the estimated SFR, the local supernova discoveries 
strongly constrain any physical effects that might increase the energy output of massive stars, including, but 
not limited to, stellar rotation. The cosmic SFR and cosmic supernova rate data on the other hand show 
the opposite trend, with the cosmic SFR almost double that inferred from cosmic supernovae, constraining a 
significant decrease in the energy output of massive stars. Together, these lines of evidence suggest that the 
true SFR calibration factors cannot be too far from their canonical values. 
Subject headings: stars: rotation - stars: formation - supernovae: general 



1. INTRODUCTION 

The star formation rate (SFR) is one of the principal pa- 
rameters characterizing galaxies and their evolution. Much 
effort has been put into deriving SFR indicators from lu- 
minos ities spanning w a velengths from the radio to the UV 
(e.g.. iKennicut^ 119981: [Be3 120031: IMoustakas et al] 1 2006 



Calzetti et a l. 2007; Salim et al. 2007: ' Kennicutt et al.ll2009 : 
Rieke et al.li2009; .Calzetti et al...2010.) . The efficacy of SFR 
indicators depends on being able to model the properties of 
stellar populations at those specific wavelengths. The model- 
ing is typically condensed into "calibration factors" that are 
the ratios of the SFR and the luminosity in a given band, for a 
given assumed underlying stellar population. For example, 
the SFR estimated using the Ha indicator is the observed 
Ha luminosity, appropriately dust-corrected, multiplied by 
the Ha calibration factor 

Calculating the calibration factors for the various SFR in- 
dicators requires a great deal of complex physics, including 
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challenging problems such as the evolution of massive stars 
and its dependence on stellar parameters such as mass, metal- 
licity, and rotation: stellar atmospheres: the binary fraction 
and the resulting interactions: the shape of the stellar ini- 
tial mass function (IMF): whether star formation is contin- 
uous or variable; and so on. The status was summarized by 
iKennicuttI ( Il998l) . who indicated a spread of ~ 0.1 dex in 
Ha calibration factors (^ 0.3 dex in the UV) in the litera- 
ture when compared for common assumptions of metallicity 
and IMlQ. These uncertainties have typically been subdomi- 
nant compared to the far larger spread between SFR measure- 
ments made by different surveys and indicators. However, the 
global SFR density and its redshift evoluti on have been in- 
creasingly well mapped in recent years (e.g., lLillv et al.|[r996l ; 
lMadau et al.|[T996l ll998l: ISteidel et alJ[T999l:lHopkinsll2004 . 
The spread between different SFR density measurements at 
low redshifts < z < 1 is n ow only several tens of per- 
cent dHopk ms & Beacoml26ol hereafter motivating 
a renewed look at the calibration uncertainties. 

One of the potentially very important physical inputs not 
considered in many early calibration studies is stellar rotation. 
Although the standard non-rotating theories of stellar evolu- 
tion have been highly successful, new data lead to a number 
of serious discrepancies, and the quantitative effects of rota- 
tion have been significantly advanced in the past decade (see. 

Reasonable variations of the IMF introduce an overall uncertainty of a 
factor of ~ 2, However, we will test the SFR using supemovae in a way that 
is weakly affected by changes to the IMF; see Section 4, 

' HB06 also assumes a common metallicity and IMF, and notes that IMF 
vaiiations dominantly affect the accuracy (normalization) instead of the pre- 
cision (scatter), as long as the IMF variations are universal. 
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e.g., reviews by iMaeder & Mevnelll2000t lLangeij|2012b . Im- 
portantly, stellar rotation impacts stellar e volution through its 



effect on stellar mixin g and ma ss loss (e.g ., Meynet & Maeden 



I2OOOI: iHeger & Lan "ger 2000: .Brottet al1l201lb . This affects 



the massive O and B stars that drive the determination of the 
SFR calibration fac tors. 

Recently, Ekstr om et al.l (1201 2b presented evolutionary 
tracks of single stars at solar metallicity, including a com- 
plete grid of tracks that have a detailed treatment of axial ro- 
tation (the "ROT" tracks, to distinguish from the non-rotating 
"NOROT" tracks). In the ROT tracks, they assume initial ro- 
tations of Vini — 0.4ucrit, whcrc Ucrit IS the velocity where the 
centrifugal and gravitational accelerations are equal, and the 
ratio of 0.4 is motivated by the peak of the veloci ty distribu- 
tion observed in young B stars (lHuangetal.ll2010l) . 

If rotation a t the level assumed in the ROT tracks of 
lEkstrom et al.l (l2012b is ubiquitous among stars, it implies 
a number of changes and improvements over the NOROT 
tracks. For example, they provide a better description of the 
observed properties of Wolf-Rayet stars, including the Wolf- 
Rayet type ratios, as well as the resulting core- collapse su- 
pernova (CCSN) type ratios (iGeorgv et al.ll2012l) . They also 
show considerably improved agreement with the positions of 
evolved yellow and red supergiants on the HR diagram, as ob- 
served in M 33 (Drout et apl2012b and the Large Magellanic 
Cloud (iNeugent et al.l2012r Similarly, it has been shown that 
the increased mass loss rates result in better agreem ent with 
the o bserved yellow supergiant CCSN proge nitors (iGeorgv 



2012 ). Furthermore, using the S tarburst99 (iLeitherer et al 



1999 ) evolutionary synthesis code. lLevesque et al.l ( 2012 ) has 



shown that the ionizing radiation field produced by a popula- 
tion of ROT stars is considerably harder than those produced 
by a population of NOROT stars. The difference reaches an 
order of magnitude in the ionizing A < 228 A flux. They also 
note that the bolometric luminosity of the ROT population is 
higher by a factor ^ 2 for the first ~ 10 Myr 

Such changes to the photon output of stellar populations 
will impact SFR measurements that are derived from observed 
galaxy luminosities. We quantitatively assess whether these 
proposed changes are consistent with observations of star for- 
mation and CCSN rates. 

We quant ify for the first t i me th e effects of stellar rotation as 
modeled by Ekstrom et al. (|2012b on SFR calibration factors. 
Pr eviously, Leithererl (|2008|) used the rotational stellar tracks 
of iMevnet & Maedeil ( I2OOOI) and reported ^ 25% reductions 
in SFR calibration f actors. We base our calculations on the 
newer ROT tracks o f lEkstrom et all (1201 2|) and find larger ef- 
fects. 

We then quantitatively investigate whether the new SFR 
calibrations are consistent with recent SFR and CCSN data. 
Since CCSNe mark the ends of the most massive stars formed 
in a star formation burst, their occurrence is an excellent proxy 
for recent star formation. CCSNe are particularly suited for 
our purposes of testing the effects of the ROT tracks on SFR 
calibrations, because CCSNe provide a benchmark SFR that 
is minimally affected by stellar rotation: in other words, while 
stellar rotation is a major parameter affecting what type of 
CCSN a collapsing star becomes, the total number of CC- 



SNe, when all Types II and Ibc are included, likely remains 
only weakly affected. In addition, the SFR-CCSN compari- 
son is only weakly dependent on the IMF, which allows us to 
avoid the large uncertainty introduced by IMF variations. 

We focus on both the cosmic and more local distance 
regimes. For the latter, we make new calculations of the SFR 
in the nearby Universe within 1 1 Mpc (we henceforth refer 
this volume simply as the "local volume"), update the dis- 
covery list of local CCSNe, and perform a comparison of the 
two. The local volume provides a unique test that is comple- 
mentary to comparisons at cosmic distances, since it is depen- 
dent on different inputs with different systematic effects. The 
completeness of CCSN searches has also improved in recent 
years, enabling a meaningful analysis with decent statistics. 
For example, CCSNe in the local v olume have been rece ntly 
utilized in the study of d im CCSNe dHoriuchi et al.ll201 lb . flie 
mass range of CCSNe dBotticella et al."|2012|), and strongly 
dust-obscured CCSNe (iMattilaet al. .201^. Here, we use it 
to study the SFR calibration factors. 

The paper is organized as follows. In Section|2l population 
synthesis results are presented and the SFR calibration factors 
are calculated. In Section[3] the effects on SFR measurements 
are discussed, for both the cosmic and local distances. We 
then study the comparison with observed CCSN rates in Sec- 
tionlD We finish with discussion and a summary in Section|5] 
Throughout, we adopt the standard ACDM cosmology with 
Vtm = 0.3, 17a = 0.7, and i?o = 70 km s"! Mpc^^. 

2. EFFECTS OF NEW STELLAR TRACKS ON SFR ESTIMATORS 

2.1. Integrated properties of stellar populations 

We use the PEGASE.2 evo lutionary synthesis code 
dFioc & Rocca-Volmerangel 1 19991) to investigate integrated 
properties of synthetic stellar populations . The default tracks 
come mainly from t he Padova group ( Alongi e t al.l 19931 
Bressan et all 1 19931 iFagotto et all [l994a.b.c: Girardi et al .1 



19961 henceforth "Padova96"). These are non-rotational stel- 



lar tracks covering a wide range of mass and metallicities. 
To obtain a complete set of stellar evolutionary tracks, PE- 
GASE.2 combines AGB, white dw arf, and unevolving stel- 
lar evolutio ns by var ious authors ("Groenewegen & de Jong] 
[T993: Alth aus & Benve nuto 1997; Chabrier & Baraffe 1997i). 
Throughout, we adopt a fixed metallicity of Z = 0.014 in or- 
der to match the metallicity of the tracks by lEkstrom et al.l 
(I2OI2I). -furthermore, we adopt the supernova model B of 
IWooslev & Weave3(ll995b . a close binary fraction of 0.05, no 
galactic winds, no treatment of dust extinction, and a mini- 
mum core-collapse supernova mass of SM©. In PEGASE.2, 
the supernova model and the minimum mass are only used 
for calculations of metal enrichment. We checked that chang- 
ing these parameters do not change the integrated popula- 
tion emission properties noticeably. Unless otherwise stated, 
we adopt a Salpeter IMF with slope dN/d M oc M" wit h 
a = -2.35 over the mass range O-l -lOOM© (ISalpetedfT955b . 

The NOROT and ROT tracks of lEkstrom et al.l(l20I2b are 
for fixed metallicity {Z = 0.014) and extend to the core He- 
flash, early-AGB, and the end of the core C burning for low, 
intermediate, and massive stars, respectively. Since we are 
interested in the first ^ 100 Myr of evolution and in bands 
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Fig. 1. — Top panel: bolometric luminosity ((iiJboi/dt = -Lbol) following 
a burst of star fonnation with a fixed stellar mass of ID^Afo, for the stellar 
tracks of Padova96 (non-rotating; dotted), GenevaOl (non-rotating; dashed), 
and for Ekstrom et al. (2012): NOROT (non-rotating, dot-dashed) and ROT 
(rotating, red solid). Bottom panel: the luminosity per logaiithmic time bin 
(dE^i^l/dlog t), illustrating the nearly constant energy output after reaching 
a peak at ^ 4 Myr when the most massive stars begin to disappear due to 
core collapse. 

where massive stars dominate, these tracks are generally suf- 
ficient. For example, the central evolution of massive stars 
beyond C-burning is usually short enough that the surface 
properties are no longer modified, and we also do not need 
to include white dwarf tracks. We adopt the same supernova, 
dust, binary, and IMF assumptions as described above for 
the Padova96 population. For com parison purposes, we also 
consider the non-rotating t racks bv ISchaller et al.l (Il992l) and 
eune &Schaered(l200ll) . henceforth "GenevaOl". 

We first simulate an instantaneous burst of star formation 
with a fixed stellar mass of lO^M© and follow the evolu- 
tion of this stellar population up to 100 Myr Such a co- 
eval stellar population can be used for studying different stel- 
lar samples as the population ages. Figure [T| shows the re- 
sultin g bolometric luminosity. We confirm previous w orks 
(e.g., iLeitherer & EkstromI 1201 it iLevesque et al.l I2OI2I) and 
find that the ROT population is 0.4-0.5 mag more luminous 
than the NOROT population. This difference is a consequence 
of rotation-induced mixing. This allows stars to bum more 
hydrogen and increase the main sequence lifetimes by about 
25%. They also have higher lu minosities than th e NOROT 
stars. These echo the results of iBrott et al.l (1201 lb who also 
show luminosities and temperatures increase for fast rotation. 
We continue to observe this difference out to 100 Myr, which 
is consistent with the fact that the rotational effects are ob- 
served for all stars above ^ 2Mq (lEkstrom et al.ll2012i) . In 
the bottom panel we show how the energy output appears in 
bins of logarithmic time. 

The NOROT population evolves similarly to those based on 
the older non-rotating tracks of Padova96 and GenevaOl, but 
after several Myr becomes sUghtly underluminous. This is 
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Fig. 2. — Quantities following a burst of star fomiation: UV (0.2/jm) lumi- 
nosity (top panel). Ha luminosity (middle panel), and the CCSN rate (bottom 
panel). Star formation burst properties are same as in Figure □ The UV lu- 
minosity evolves similarly to the bolometric luminosity, since a wide range 
of stars (~ 3-2OM0) contribute to the UV band. On the other hand, the Ha 
luminosity is powered by the most massive stars (> 2OA/0) and its behavior 
is qualitatively different. 

due to the updated overshooting and mass-loss prescriptions 
of the NOROT tracks. These decrease the main sequence life- 
time of M < 3OA/0 mass stars and increases the main se- 
que nce lifetime of M > 3 0Mp, stars, respectively (see Figure 
6 of lEkstrom etaLll2012l) . Therefore, after a few Myr when 
the most massive stars have exploded, the difference between 
populations are dominated by M < 30Mq stars, and the 
NOROT population is expected to become dimmer compared 
to the Padova96 and GenevaOl populations. However, the dif- 
ferences are small, and we will henceforth quote effects of 
rotation relative to the updated non-rotational NOROT tracks 
unless specifically stated otherwise. 

In Figure 121 we show the UV (0.2/im) luminosity, the Ha 
luminosity, and the CCSN rate following the same star for- 
mation burst. The UV luminosity evolution is similar to that 
of the bolometric luminosity, because the dominant contribu- 
tion comes from a wide mass range of stars; for our adopted 
tracks, ^ 3-20Mq stars. On the other hand, the Ha luminos- 
ity shows a qualitatively different behavior This is because 
the Ha emission is largely sensitive to photons shortward of 
the Lyman limit that dominantly arise from massive young 
stars with M > 20 Mq. The Ha thus falls when these massive 
stars disappear in core collapse. The progenitors of CCSN are 
M > 8Mq stars that fall in between the UV and Ha in terms 
of mass range. Thus it rises when the most massive IOOMq 
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Fig. 3. — Same as Figure |2]but for a continuous stai' formation rate of 
50 Mq yr^^. The Ha reaches equilibrium earliest since they are sourced by 
the most massive stars only. The progenitors of CCSNe are M > SMq stars 
that reach equilibrium later, while the UV has some contribution from stars 
with lifetimes longer than 100 Myr. 

Stars core collapse and falls when the least massive 8Mq stars 
undergo core collapse. For the CCSN rate, the difference be- 
tween tracks is only due to the different main sequence life- 
times of the CCSN progenitors. For example, the ROT tracks 
with their longer main sequence lifetimes shows a later rise in 
the CCSN rate. On the other hand, the UV and Ha properties 
also strongly reflect the differences in luminosity and surface 
temperature. 

2.2. Star formation rate calibration factors 

A variety of indicators are used as tracers of recent star 
formation activity, all directly or indir ectly probing the cur- 
rent 



j opulation of massive stars (e.g.. | Kennicutt 



2003 



2007 



Moustakas et al.ll2006l: ICalzetti et al.ll2007l: 



19981 iBeii 



SaUm et al 



kennicutt et aiT (2009l: iRieke et alj|2009l : ICalzetti eTal 



2010|). The UV stellar continuum and the Ha nebular emis- 
sion have traditionally been extensively used due to their ac- 
cessibility. We compute the UV and Ha calibration factors 
using PEGASE.2 assuming standard assumptions of constant 
star formation rate for at least 100 Myr This is the continuous 
star formation approximation, which provides enough time to 
allow the birth and death of massive stars that dominate the 
luminosity in the UV and Ha to equilibrate. It is applica- 
ble at least at low-z where the fraction of the total SFR con- 
tributed from young starbursts is thought to be < 10% (e.g., 
ISalim et al...2005,) . 
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Fig. 4. — UV (top panel) and Ha (middle panel) calibration factors cal- 
culated under the continuous star formation approximation stalling at time 
t = 0. Calibration factors are defined such that SFR = f^L X, I.e., in 
units of (M0yr~^)/(ergs~^ Hz~^) and (M0 yr~-'^)/(ergs~-'^) for the 
UV and Ha, re spectively. For comparison, we mark the values quoted in 
iKennicutil i 19980 as a data point. In the bottom panel we show the CCSN 
calibration factor (see Section|4]for details). 

In Figure [3] we show the UV (0.2/im) luminosity, the Ha 
luminosity, and the CCSN rate, all for a continuous star for- 
mation rate of 5OM0yr~^. As can be seen, the Ha lumi- 
nosity and CCSN rate comfortably reach equilibrium before 
100 Myr By contrast, the UV is still marginally increasing 
because the UV has contributions from stars with main se- 
quence lifetimes longer than 100 Myr However, this is not a 
large fraction. The UV luminosity only increases by another 
~ 15% over the next few hundred Myrs before reaching equi- 
librium. It should be mentioned that the CCSN rate for differ- 
ent stellar tracks all tend to the same value. This is because 
we have defined the CCSN rate by the mass range of stars 
undergoing core collapse (see Section|4]for more details). 

We define the calibration factor by 



SFR: 



(1) 



where SFR is in units M0 yr 

'1 



fxLx, 

^, and Lx is the dust-corrected 
luminosity in units crgs^^ for the Ha and ergs^^ Hz^^ for 
the UV. Results for the Padova96, GenevaOl, NOROT, and 
ROT tracks are shown in Figure |4] where their time evolu- 
tion are shown, and their values for 100 Myr after the onset 
of star formation are summarized in Table [U Firstly, we see 
that the new N OROT cahbrations are similar to the Padova96, 
GenevaOl, and IKennicutil l fl99i) calibrations, although they 
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TABLE 1 

Summary of Z = 0.014 calibration 

FACTORS FOR SALPETER IMF AT 100 MYR 
(CONTINOUS SFR) 
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NOROT (2012) 
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ROT (2012) 
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Note. — Calibrations defined as /^^ = 
SFR/Lj^. The Ha calibration has units of 
{Mq yr~-'^)/(erg s"-'^) and the UV calibration has 
units of (A-/Qyr-l)/{ergs"l Hz-^). 



are slightly lar ger For example , compared to the commonly 
used values of lKennicuttl(ll998l) . the UV calibration is larger 
by 10%. Secondly, the ROT tracks yield smaller calibra- 
tions. For example, the ROT UV and Ha calibrations are 
smaller by 30% and 40% when compared to the NOROT cali- 
brations. These follow from the increased main sequence life- 
times and lumi nosities of rotatin g stars. The widely adopted 
calibrations of iKennicutll d 19981) are shown for comparison 
with error bars th at reflect the full range of values reported by 
iKennicutl ( Il998l) . Most of our calibrations from Table [T| are 
within this range, but the ROT Ha calibration is significantly 
lower 

We change the PEGASE.2 input parameters to test the ro- 
bustness of these results. One of the most important inputs 
for calculating the calibration factor is the IMF. Generally, a 
shallower (steeper) IMF results in relatively more (less) mas- 
sive stars and smaller (larger) calibration factors. We consider 
three IMFs for comparison: the "SalA" IMF, as proposed by 
iBaldrv & Glazebrookl (l2003h . with the Salpeter shape above 
0.5Mq (high-mass gradient of —2.35) and a suppression be- 
lo w 0.5M(7i (low-mass gr a dient o f —1.5); the "BG" IMF, also 
bv iBaldrv & Glazebrookl (|2003|) . with a high-mass gradient 
of -2.1 5 and a l ow-rn ass gradient of —1.5; and finally, the 
IMF of lKroupal (l200lh . The SalA, BG, and Kroupa IMFs 
respectively yield UV calibration factors that are ^ 0.77, 

0.50, and 0.83 times those of the Salpeter IMF, afl for the 
NOROT tracks. The Ha changes slightly more: by ~ 0.77, 
~ 0.42, and ~ 0.79, respectively. 

Repeating the exercise with the ROT tracks, very simi- 
lar changes are observed. For example, the SalA, BG, and 
Kroupa IMFs yield UV calibration factors that are ^ 0.78, 

0.49, and ^ 0.83 times those of the Salpeter IMF, and Ha 
calibration factors that are ^ 0.77, ~ 0.41, and ^ 0.79 times 
those of the Salpeter IMF. Therefore, while the absolute val- 
ues of the calibrations change as a result of the IMF shape, the 
relative importance of stellar rotation is always a reduction of 
30-40%. We also stress here that the IMF does not strongly 
impact the comparison with the CCSN rate (see Section^. 

Changing the mass range of the IMF also makes a differ- 
ence to the SFR calibrations. For example, increasing the 
high-mass cutoff from our canonical IOOMq to 120Mo leads 
to more massive stars (for a fixed total stellar mass) and a 3% 
and 15% decrease in the UV and Ha NOROT calibrations, re- 



spectively. Some authors adopt an upper limit of 60 Mq which 
would have the opposite effect: 12% and 60% increases, re- 
spectively. We note here that the use of 6OM0 necessarily 
underestimates the Ha calibration because the important mas- 
sive stars are not included. Once again, the effect of rotation 
is a 30-40% reduction in the calibration factors regardless of 
the mass range. This is partly due to the assumption of a rota- 
tion velocity i)i„i = 0.4wcrit for all stellar masses of the ROT 
tracks. 

Metallicity is also important since the evolution of massive 
stars are strongly aff ected by metal-driven mass loss (e.g., 
iKudritzki et al.l [l989l) . Generally, low-metallicity stars are 
larger and more luminous, leading to smaller SFR calibration 
factors. For the Padova96 tracks and a Salpeter IMF, changing 
the metaflicity to Z = 0.001 {Z = 0.1) causes the UV calibra- 
tion to change by —2% (+10%) when compared to our canon- 
ical Z = 0.014 estimates, and the Ha calibration to change 
by —30% (+50%). We do not quantitati vely study rota t ion ef- 
fects for different metallicities since the lEkstrom et al.l ( |20I2|) 
tracks are only for a single metallicity. However, rotation ef- 
fects are less important for high-metallicity stars whose evo- 
lutions are more s trongly affected by metal-driven mass loss 
dBrott et al.ll201 lb . Other parameters such as the minimum 
mass for CCSNe, the supernova chemical abundance model, 
and the presence of Galactic winds, only slowly affect the lu- 
minous outputs of stars through feedback of chemically en- 
riched gas. 

Other star formation indicators that are routinely used in- 
clude the [Oil] A3727 forbidden line, infrared, and radio emis- 
sions. We do not calculate calibrations for these directly. 
However, because they are often associated with either the U V 
or Ha, we can discuss how they are affected by stellar rota- 
tion. For example, the [Oil] luminosity is not directly coupled 
to the ionizing luminosity, but as a SFR tracer it is empirically 
cali brated based on Ha, often using samples of nearby galax- 
ies (|Kennicut3 [1991; [MouiSil^^ This is non- 
trivial since the [Oil] to Ha ratio depends on the luminosity, 
metallicity, and also obscuration. But it implies that changes 
to the Ha calibration factor will be carried through to a change 
in the [Oil] calibration factor 

The infrared luminosity is an indirect tracer of star for- 
mation arising from UV continuum that is absorbed and re- 
radiated by dust. It is complementary to the UV and Ha that 
that are absorbed by dus t, and calibrations have been inves- 
tigated by m any authors (Sauvage & Thuan"'l992VBell"2003j; 
ICalzetti et al. 2007; Rieke et al. 2009 ; Calzetti et al. 2010). It 
requires sufficient optical depth of the dust in the star forma- 
tion region, such that the Lir provides a calorimetric mea- 
sure of the star formation. However, in practice, the in- 
frared is contaminated by dust heated by old stellar popula- 
tions. For spiral galaxies the contaminatio n to the far-infrared 
(1-10 00/^tm) can be as high as 70% (e.g.. lLonsdale & Heloul 
Il987h . and modern calibrations therefore utilize shorter wave- 
lengths where contamination is less of an issue dCalzetti et alj 
I2OIOI) . So, rotation effects on the low-mass stellar popula- 
tion must also be considered to some degree, depending on 
the galaxy sample an d exact wavelength used. According to 
lEkstrom et al.l (|2012|) . rotation affects low mass stars (0.8- 
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2Mq) less than the massive stars, which means that the in- 
frared calibrations will be less affected by stellar rotation than 
the UV or Ha. 

Radio continuum emission is another indirect tracer of star 
formation. It is insensitive to dust attenuation, but it relies 
on the complex and poorly understood physics of cosmic ray 
production and confine ment i n super nova remnants. The 1.4 
GHz SFR caUbration of Hsl (l2003b partly avoids this diffi- 
culty by empirically using the observed tight radio-FlR corre- 
lat ion, and ca l ibratin g to be consistent with the FIR calibration 
of lKennicut3(ll998l) . This calibration is a factor ~ 2 smaller 
than the calibration o f lCondonl(ll992h that is calculated from 
first princip les. Rotation will reduce the radio calibrations of 
iBelll (12003 ) through the infrared calibration, but will not affect 
those of ICondo^ jl992l) . 

To conclude, recent modeling of stellar rotation imply re- 
ductions in SFR calibration factors by 30% for the UV and 
40% for the Ha, where we have compared the ROT popula- 
tion to the NOROT population. Comparing the R OT calibra- 
tions to the widely-used non-rotating calibrations of lKennicuttI 
[\99&) . the reductions are 10% and 40% for the UV and Ha, 
respectively. 

3. IMPLICATIONS FOR SFR ESTIMATES 

The reductions to the SFR calibration factors due to stel- 
lar rotation are somewhat smaller than the reductions caused 
by variations in the IMF 0.77 for SalA and - 0.50 for 
BG). Nevertheless, they have testable consequences. As we 
discuss in Section]?] observations of CCSNe provide a bench- 
mark SFR whose comparison with the directly measured SFR 
is minimally dependent on the IMF. Also, because stellar ro- 
tation affects the Ha indicator more than the UV, we can use 
the UV/Ha ratio as a test. 

3.1. SFR at cosmic distances 

At cosmic distances, the SFR has been measured by many 
groups using a variety of indicators. Comparisons of data 
are complicated by different galaxy selection biases, differ- 
ent assumptions of dust attenuation, diff erent luminosity cuts , 
cosmic variance, and other issues (e.g., iHopkins et al.ll200ll: 
[Hopkin s 2004). To address these, HB06 selected recent SFR 
data, converted them to a common cosmology, dust correc- 
tion scheme, IMF, and calibration assumptions. Using their 
compilation, they found a best-fit cosmic SFR density with a 



z = normalization of 2.0toi ^ 10"^ Mpc-^ ^yr"^ (la 



errors) for a Salpeter IMF and SFR calibrations of lKennicuttI 
(Il998l) . 

The HB06 best-fit uncertainty is ±20-30% in the most 
constrained redshift range of < z < 1. In this redshift 
range, the HB06 fit is based mostly on UV-derived SFR mea- 
surements, except f or one Ha measurement at z w 0.01 by 
iHanish et al. I (l2006h . Cha nging the UV calibration from their 
adopted Kennicutfl ( 1 1 9981) values to our ROT values only re- 
sults in a mild 10% reduction in the SFR, which is smaller 
than the fit uncertainty. At larger redshifts, a wider range of 
indicators are used in the HB06 fit. Since stellar rotation af- 
fects different indicators by different degrees, stellar rotation 
will contribute additional spread as well as reduce the SFR 



normalization. However, the HB06 best-fit uncertainty is fac- 
tors of 2 or more at these larger redshifts, so rotation would 
not be a dominant contributor to the uncertainty. 

3.2. SFR at local distances 

In contrast to the cosmic SFR, studies of the local vol- 
ume within ~ 10 Mpc have the advantage that all use almost 
identical galaxy samples. The 11 Mpc Ha and Ultraviolet 
Galaxy Survey (IIHUGS) is designed to provide a census of 
the SFR in the volume within 1 1 Mpc using the Ha indicator 
jKennicutt et alj |2008|) . The distance limit was set to have 
both statistically significant galaxy numbers (436 galaxies) 
and galaxy completeness (Mb ~ —15 mag at the distance 
limit 11 Mpc; however, the number of galaxies as a function 
of distance does not rise as fast as expected from a constant 
density of galaxies, as shown in Figure|5]). These local galax- 
ies have also been observed by the GALEX and Spitzer 
satellites, providing multi-band UV and IR photometry: these 
are the GALEX Nearby Gal axy Survey (390 galaxies, includ 
ing 3 63 from the 1 IHUGS: loil de Paz et al.ll2007l: IE 



.^ee et 



iD 



I20T U) and the composite Local Volume Legacy survey (LVL, 
with a total of 25 8 galaxies, including 231 from the 1 IHUGS; 
iDale et al.ll2009l) . 

Using a sample of 300 local galaxies, iLee et alJ (l2009l) 
carefully investigated the Ha- and UV-SFR. After correct- 
ing for dust, they find that the Ha-SFR tend to be smaller 
than the UV-SFR, log(SFRuv/SFRHa) ~ 0.13, i.e., a fac- 
tor of 1.35, for the most star-forming galaxies with Ha-SFR 
> lO-i '^Mgyr-^ In another studv. iBotticella et al.l (1201 2|) 
studied 312 local galaxies and derived dust-corrected Ha and 
UV luminosities of Lhq = (99 ± 5) x lO^^ergs"^ and 
Luv = (88 lb 6) X lO^^^ergs"^ Hz~^, respectively. Using the 
calibrations of lKennicuttI jl998h . they find Ha- and UV-SFR 
values of 78 ± 4 Moyr"^ and 123 ± 8 M0yr~^, respectively. 

If we appl y the ROT calibration factors to the local lu- 
minosities of Botticella et al. I (l2012l) . the difference between 



the Ha- and UV-SFR will widen from their nominal ra- 
tio of UV-SFR/Ha-SFR « 1.6 to UV-SFR/Ha-SFR > 
2. This difference is considerably larger than thes reported 
measurement errors, and disfavors ubiquitous stellar rota- 
tion in the local galaxies. However, it is premature to con- 
clude that stellar rotation is not allowed by local SFR es- 
timates. On an individual galaxy basis, different indica- 
tors give dis crepant SFR estimat es, sometimes up to fac- 
tors of 10 dHopkins et al.ll2003b . Although the difference 
is much reduced for larger samples of galaxies, systematic 
uncertainties in flux measurements, dust attenuation correc- 
tions, stochasticity in the formation of high mass stars, vari- 
ations in the IMF, and departures from standard recombi- 



1 — 1 

tween the two indicators (e.g. Buat et al. 


1987; 'Buat 


1992; 


Hopkins et al.ll200ll; Buat et alji2002: Bell & Kennicutt 2001: 


Islesias-Paramo et al.ll2004l; Sahm et al.ll2007: 


Bothwell et al. 


2OO9I). Of these, the most plausible source of problems is dust 



discrepancy arises due to uncertain dust corrections of atypi- 
cal galaxies that dominate the SFR budget; for example, high 
inclination galaxies (e.g., M 82) or galaxies with large fore- 
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ground Galactic extinction (e.g., NGC 6946). 

3.3. Total SFR within 11 Mpc 

We derive the total SFR within the 11 Mpc volume and 
compare to the cosmic SFR density. For this we must as- 
sume a lo cal overdensity fac t or Fo r example, at distances of 
1-5 Mpc, iDrozdovskv et al.l (l2008h demonstrates more than 
factor 2 overdensity over the cosmic SFR. iKarachentsev et afl 
( I2OO4I) studied the B-band luminosity density of the local 8 
Mpc and found it to be 1.7-2.0 times the global luminosity 
density as measured by the SDSS and the Millennium Galaxy 
Catalogue. Using an overdensity factor of 1.7, and adjust- 
ing to the NOROT calibration, the z ~ extrapolation of 
the HB06 best-fit cosmic SFR density predicts a total SFR 
within 11 Mpc of 220^30 Moyr'^ Similarly, the ROT cali- 
bration yields 170^3q M0yr~^. In both estimates we use the 
UV calibration, since the HB06 compilation in the redshift 
range < z < I is dominantly based on the UV indicator. 

We directly estimate the total SFR within 11 Mpc by start- 
ing with the IIHUGS + GALEX galaxy catalog and apply- 
ing the same two-tier Galactic latitude cut as is used in the 
LVL survey, namely, |6| > 20° for distances < 3.5 Mpc and 
|6| > 30° for distances of 3.5 to 11 Mpc. This leaves a sam- 
ple of 282 galaxies with positive Ha and FUV flux measure- 
ments. The fraction of the 11 Mpc sphere included after these 
cuts is ~ 51%, and we estimate the total SFR by summing the 
SFR of the galaxies and dividing by the volume fraction. 

The Ha luminosities are corrected for the [Nil] line con- 
tamination, underlying stell ar absorption, and Gala ctic fore- 
ground extinction following iKennicutt et al.l (|2008|) . We cor- 
rect for internal attenuation using the empiric al scaling correc- 
tion with the host galaxy i?-band magnitude dLee et all2009h 
of 



A 



Ho 



0.10 Mb > -14.5 

1.971 + 0.323Mb + 0.0134Af| Mb < -14.5 

(2) 

For the FUV measurements, we correct for internal attenu- 
ation based on the tigh t correlation with the TIR/FUV ratio 
(iBurgarella et al ] l2005h . 

A(FUV) = -OmSx^ + 0.392x^ + 1.094a; + 0.546, (3) 

applicable for galaxies with recent star formation such as the 
spirals and irregulars. Here, x = log[L(TIR)/L(FUV)obs]- 
When the TIR data gives a negative correction, or when TIR 
data is not ava ilable, we u se the empirical relation ^(FUV) — 
1.8A(Ha) of lLee et all (120091) . These corrections are sim- 
ilar to thos e applied in previous studies (iLee et al. 1 12OO9I: 
iBotticellaet al. 2012J. Finally, we modify the UV calibrations 
of Table [l] to reflect the central wavelength of the GALEX 
FUV filter, 0.1532^m. For the NOROT, we find that 



logSFR[Moyr-i]=log(L 



FUVxorr 



[L0])-9.51, (4) 



which is very similar t o the G ALEX calibration reported in 
esias-Paramo et al.l (l200q) who uses the Starburst99 syn- 



thesis code (iLeitherer et al 



I999I) with the solar metallicity 



TABLE 2 

Summary of total SFR estimates within 11 
Mpc 



SFR (Mq yr-i) 


Cosmic extrapolation (NOROT) 


220 ± 40 


Direct Ha (NOROT) 


120 ± 20 


Direct UV (NOROT) 


170 ± 20 


Cosmic extrapolation (ROT) 


170 ± 30 


Direct Ha (ROT) 


80 ± 10 


Direct UV (ROT) 


130 ± 20 


CCSN witli 08S-like objects 


360+;,f 


CCSN without 08S-like objects 


27Q+110 



Geneva non-rotating tracks and a Salpeter IMF from 0.1- 
1OOM0. 



Note. — These are the total SFR, i.e., sky-coverage 
con'ected where appropiiate (see text). All for a 
Salpeter IMF. The cosmic extrapolation is corrected for 
the local overdensity by a factor of 1.7. For CCSNe, 
a supernova mass range of 8-40 Mq is adopted; the 
upper mass limit is not important (see text). 



We summarize our results in Table |2] The direct SFR es- 
timates are generally smaller than the cosmic extrapolations, 
but the UV-SFR is consistent within the uncertainties. The Ha 
estimates, on the other hand, fall short of the cosmic extrapo- 
lation. One possibility is that we have slightly overestimated 
the local overdensity, since we adopt an overdensity derived 
from ^ 8 Mpc scales for the entire 11 Mpc volume. 

As mentioned in the previous section, the stochastic nature 
of star formation can be an issue. Since PEGASE.2 is a deter- 
ministic synthesis code, it does not cover the stochastic behav- 
ior at low SFR rates. This occurs when the number of massive 
stars falls below ^10. However, we are summing over a suf- 
ficient number of galaxies with a sufficiently large total SFR 
that the stochasticity is not a problem. More quantitatively, 
the Ha indicator is related to stars with M > 20 Mq. Adopt- 
ing the Salpeter IMF, 5 x 10'^ ALq of stellar material must be 
made for there to be 10 such stars, and this must continue 
at least over the lifetime of the massive stars, some 3 Myr 
Thus, the required SFR is at least 2 x 10"'^ Mq yr~^ per 
galaxy. If we adopt more modern IMFs where the number of 
low mass stars are suppressed relative the Salpeter IMF, there 
will be more massive stars for a given total stellar mass and 
the required SFR will only decrease. 

4. COMPARISON TO CORE-COLLAPSE SUPERNOVA RATES 

The most massive stars in a star formation burst evolve the 
fastest and then explode as CCSNe. Therefore, observations 
of CCSNe are excellent proxies of recent star formation that 
is far less sensitive to modest dust. Different types of CC- 
SNe probe the SFR on different time scales; treated together, 
CCSNe probe the SFR on time scales similar to the UV indi- 
cator. A well-measured CCSN rate yields an estimate of the 
SFR normalization, and even the observation of a few CCSNe 
yields lower limits. We use both the cosmic CCSNe and CC- 
SNe in the local volume and discuss how compatible they are 
with the rotation calibrations. 

The CCSN rate, A^ccsn, is related to the SFR, A/*, by 
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TABLE 3 

Local CCSNe (in bold, total 17) and possible CCSNe (in italic, total 5) within 11 Mpc During 2000 to 2011 

Inclusive 



SN 


Galaxy 


Type 


D (Mpc 


R.A. 


Decl. 


b (deg) 


11 HUGS 
Full LVLcut 


LVL 


SN zOOzap 


NGC 0628 


IcPec 


7.3 


01 36 41.7 


15 46 59 


—45.7 


Y 


Y 


Y 


SN 2002bu 


NGC 4242 


SN?(i) 


7.4 


12 17 30.1 


45 37 08 


70.3 


V 


V 


Y 


SN 2002hh 


NGC 6946 


IIP 


5.9 


20 34 52.3 


60 09 14 


11.7 


V 
I 




" 


aJN zUUJgd 


JNuC [)o2q 


IIP 


7.3 


01 3d 41.7 


15 46 59 


—45.7 


Y 


Y 


Y 


aJN 2UU4am 


JNCjC jUJ4 (M 82) 


IIP 


3.5 


09 55 52.2 


Cf\ A f\ A'7 

69 40 47 


40.6 


Y 


Y 


Y 


SIN 2UU4aj 


JNLrC 2403 


IIP 


3.2 


r\'7 c: 1 A 

U7 -3D 51.4 


oc r\r\ 

65 36 09 


29.2 


Y 


Y 


Y 


aJM 2UU4et 


JNLrC oy4D 


IIP 


5.9 


on OA c;o q 
ZU v>4 02.v> 


60 09 14 


11.7 


Y 










TTT3 
llr 


J.O 


1 Q nc; 07 
io UO Z ( .D 


/lo nfi 
— 4y Zo UD 


io.o 


Y 


. 








Ic 


Q A 




a'i t:i or: 
— UO OL ZO 


26 1 


Y 










IIP 


8.0 


2Q ^2 7 

-LO iii7 ^.JZJ . ( 


47 11 43 


68.6 


Y 


Y 


Y 


SN 2007gr 


NGC 1058 


Ic 


9.2 


02 43 29.9 


37 20 27 


-20.4 


Y 






5^20055 


NGC 6946 


SN?(2) 


5.9 


20 34 52.3 


60 09 14 


11.7 


Y 






SN 2008ax 


NGC 4490 


lib 


8.0 


12 30 36.1 


41 38 34 


74.9 


Y 


Y 


Y 


SN 2008bk 


NGC 7793 


IIP 


3.9 


23 57 49.7 


-32 35 30 


-77.2 


Y 


Y 


Y 


NGC300-OT 


NGC 0300 


SN?(3) 


2.0 


00 54 53.5 


-37 41 00 


-79.4 


Y 


Y 




SN 2008iz 


NGC 3034 (M 82) 


II 


3.5 


09 55 52.2 


69 40 47 


40.6 


Y 


Y 


Y 


SN 2008jb 


ESO 302-14 


IIP 


9.6 


03 51 40.9 


-38 27 08 


-50.9 


Y 


Y 




SN 2009hd 


NGC 3627 (M 66) 


IIP 


10.1 


11 20 15.0 


12 59 30 


64.4 


Y 


Y 


Y 


SNlOlOda 


NGC 0300 




2.0 


00 54 53.5 


-37 41 00 


-79.4 


Y 


Y 




SN 2011dh 


NGC 5194 (M51) 


lib 


8.0 


13 29 52.7 


47 11 43 


68.6 


Y 


Y 


Y 


PSNJ12304185+4137498 


NGC 4490 


SN?(^' 


8.0 


12 30 36.1 


41 38 34 


74.9 


Y 


Y 


Y 


SN 2011ja 


NGC 4945 


IIP 


3.6 


13 05 27.5 


-49 28 06 


13.3 


Y 







_N^OTE. (l ^Thompsonetal.1 )2009h:ISmith et al.l <2011|):ISzczvgiel et alj <2012bl). (^jlPrieto e t al 



Botticella et al. (2009); Pumo et al. 12009); Prieto et al. (2010); Szczygiel et al j ]2012j) . (3)Bondetal. (2009); Bergeretal. (200SJ) 
Prietoetal (2009); Kashi et al |2010); Kochanek et al. i ,201Z) , ('') ,Elias-Rosa et alj ^2010) . 



CCSN 



iMagillet alj<201lh 

= M,//ccsN, where /ccsn is calculated as 



(2008, 2009); Smith et_al 
009); 

(^) .Coitini et al. C201 1) ;. Eraser et al. 



CCSN 



/o °" Mi>{M)dM 



nome na remains uncertain (e.g.. lFrverlll999t iKochanek et al.l 
l2008h . We adopt A/,nax = 40 M© as our canonical value, 



which together with M„ 



± 1 Mq yields /ccsn 



where O.IM© and lOOM© in the numerator is the range of 
the IMF, and Mmin to Mmax in the denominator is the mass 
range of stars that explode as CCSNe. Note that Eq. (|5]l 
is the inverse of the CCSN efficiency. Due to the steep 
IMF slope, the majority of CCSNe occur near the mini- 
mum mass threshold, making Mmin the most important pa- 
rameter Its value has been statistically determined by com- 
bining 20 Type IIP su pernova progenitor observations to be 
Mq ismartt et al.ll2009b . This is consistent 



M. 



UP 



i.5 



+1 



-i.s-' 



with the highe st masses e stimated for white dwarf progeni- 
tors, 7 M© dWilliams e t al. 2009). Thus, two different ap- 
proaches seem to be converging to i\/„iin « 8 ± 1 Mq. As we 
show below, this uncertainty affects /ccsn by ^ 20%. 

The value of Mmax is more uncertain than Afniin, but is 
less important. The maximum mass from Type IIP super- 
nova progenitors i s found to be M}^^^ = 16.5 ± 1.5 A/q 
dSmartt et al.l2009b . but consideration of other types of SNe — 
including Types Iln and Ibc — yields larger values. For exam- 
ple, it is widely expected that Type Ibc SNe originate from 
evolved massive stars that have shed their envelopes, whose 
initial masses are > 25 Mq. The spatial distribution of CC- 
SNe supports that Type Ibc progenitors must be more massiv e 
than those of Type IIP SNe jAnderson & Jamesll2008l [20091) . 
Theoretically, stars with masses above Afmax ^ 40 Mq 
may promptly form black holes whose optical transient phe- 



147^29/^^0. ^ ~ 20% uncertainty. On the other hand, vary- 
ing Afniax between 25AfQ and lO OA/© only affects .fc csN 
by ^ 5%. We refer the reader to iHoriuchi et al.l (|20I ll ) for 
detailed discussions of the uncertainties in the SFR-to-CCSN 
rate connection. 

One of the advantages of CCSNe is that the SFR normal- 
ization can be probed without significant dependence on the 
IMF. This is because the mass range of stars exploding as 
CCSNe is similar to the mass range of stars powering the 
radiation used to measure the SFR. For example, we can 
write explicitly the SFR estimated from the UV and CC- 
SNe as fvvLuv and /ccsnA^ccsn- Although changes to 
the IMF affects /uv and /ccsn individually by factors of 
^ 2, their ratio remains fairly constant. For the SalA and 
BGIMFs, /uv//cc.sN = 1.09 x lO^^Oyr-i erg^i sHz and 



1.01 X 10 



-30 



'CCSN 



s Hz, respectively, which are within 



10% of the SalpeterlMF value 1.09 x lO^^o yj.-i erg-i sHz. 



4.1. Comparison at cosmic distances 

Measurements of the CCSN rate at cosmic distances have 
rapidly incre ased in the past decad e . Taking the mea- 



surements of iCa ppellaro et al.' ('1999"); 'Dahlen et al.' ("2004^; 
ICappellaro et all (2005 ); Botticella et al. (2009); Bazi n et al] 
(l2009l) ; lLi et'IiTdioi lb . the observed CCSN rate is well-fit by 
a strong redshift evolution of (1 + z)^-^ between < 2 < 1 
with a z = normaUzation of 0.72~''n T x 



10"^Mpc"3yr-i. 
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By comparison, the CCSN rate predicted from the SFR using 
the NOROT calibration and Equation ^ has a z = normal- 



ization that is twice as large, 1-51q'2 ^ 10^** Mpc^"* yr 
Uncertainties associated with calculating the CCSN rate from 
the SFR se em not to be able to satisfactorily explain the diff er- 
ence (e.g.. lHopkins & Beacomll2006HHoriuchi et al.ll20Tll) . 

Rotation helps bridge the normahzation difference by de- 
creasing the SFR and hence the predicted CCSN rate. How- 
ever, the reduction is expected to be at most ^ 30%, which 
is not enough to bridge the difference completely. Further- 
more, it has been argued that a combination of intrinsically 
dim and heavily dust-obscured CCSNe being missed in CCSN 
surveys could be res ponsible for the observed discrepancy 
jHoriuchi et al.ll201 lb. Indeed, the most recent direct CCSN 
rate measurement bv lDahlen et"aLl (|2012|) takes into account 
missed dim CCSNe and find the discrepancy with the SFR is 
substantially reduced, supporting this explanation. 

4.2. Comparison at local distances 

Discoveries of CCSNe in the nearby volume have also im- 
proved recently. In Table [3] we show a list of CCSNe dis- 
covered within ^ 11 Mpc in the past 12 years from 2000 
to 2011 inclu si ve. Si milar compilati o ns ca n be found in 
iHoriuchi et all (1201 lb : iBotticella et"!!] ( l2012h : iMattila et all 
(l2012h . We exclude p revious years due to known supernova 
incompleteness issues (iHoriuchi & Beacom"201QV Distances 
are taken from the IIHUGS catalog (Kennicuttet al. 2008). 
These are estimated from a combination of direct stellar dis- 
tances (where available) and flow-con-ected radial velocities. 
The L ocal Group flow model of Karachen tsev & Makarovl 
(Il996h is used, which the authors find to be more accurate 
for the 11 Mpc volume than a Virgocentric flow model. The 
last three columns show whether the host galaxy is included 
in the respective galaxy catalogs as labeled. 

Included in Table [3] are SN 2008S-Uke events (SN 2008S, 
SN 2002bu, NGC300-OT, SN 2010da, and PSN 
J12304185H-4137498) whose true nature continues to be 
debated. They are characterized by explosions that are 2-3 
mag dimmer than regular CCSNe, with narrow emission 
lines, and evidence of internal extinction; their progeni- 
tors are also dust enshrouded (iPrieto et all |2008|) . While 
they are relatively common a mong nearby CCS Ne, their 
progenitors are extrern ely rare (iKhan et al.l l2010l) . leading 
[Thompson et al.l (l2009b to propose that many massive stars 
go through this dust-obscured phase shortly (~ lO'' yr) before 
the explosion. SN 2008S-like events may be true CCSNe, 
perhaps of the electr on-capture type arising f rom asymptotic 
giant branch stars (|Thompson et al.l 120091 : IBotticella et al.l 
120091: iPumo et al.ll2009l) . orthev may be extrern e versions of 



non-e x plosive outbursts ([Thompson et al 
20091 BondetalJ 12009^ iBerger et al.1 



2010i:.SmithetaI 



20091: 
2009t 



Smith et al 



Kashi et al 



201 lb . For a com prehensive discu s sion o f 



these and re la ted rn od els, see, e.g .. 
ISmiflietal.1 (1201 lb : iKochanelJ 
(I2OI2I) . 



Thompson etaP (l2009b 



201 lb : iKochanek et alj 



The compilation is a lower limit of the true CCSNe occur- 
rence. This is because the CCSN discoveries are mainly made 
by surveys and amateurs with varying search strategies and 



methods, instead of a complete systematic survey. Although 
modern telescopes using CCDs are capable of discovering a 
typical CCSN within 11 Mpc, CCSNe can be missed due to 
a combination of incomplete galaxy coverage, low survey ca- 
dence, heavily dust-obscured or dim CCSNe, or obstruction 
by the Galactic plane and the Sun. Indeed, the majority of CC- 
SNe have occurred in large, well-known galaxies: all 22 CC- 
SNe occurred in the full IIHUGS galaxy catalog. Also, 15 of 
these occurred in our galaxy sample with LVL cuts. Since the 
LVL sky-coverage is only 51%, we are likely missing some 
CCSNe occurring in the remaining volume. Furthermore, the 
northern hemisphere is more closely observed, resulting in 
somewhat more CCSNe discovered with positive declinations 
than with negative declinations (14 versus 8). As a specific 
example illustrating these points, SN 2008jb remained undis- 
covered by all galaxy-targeted supernova searches despite its 
close distance of 9.6 Mpc. The host galaxy was a dwarf 
galaxy in the southern hemisphere, and the CCS N was only 
discovered in archival images of all-sky surveys (iPrieto et al.l 
l2012h . 

There are 17 definite CCSNe and 5 SN 2008S-like events 
in Table [3] With our LVL sky-coverage cuts, this is re- 
duced to 11 definite CCSNe and 4 SN 2008S-like events. 
From Poisson statistics, the la lower and upper limits of the 
11 definite CCSNe ai-e 7.7 and 15.4, which gives a CCSN 
rate of l.St'^'lyr^^ within the 11 Mpc volume, where we 
have applied a sky-coverage correction analogous to the lo- 
cal SFR calculations. The SFR required is simply the CCSN 
rate multiplied by /ccsn- For our canonical /ccsn, this 
gives 270igJ° M0 yr~^. Repeating the exercise including 
SN 2008S-like events gives a CCSN rate of 2.5toi Y^^^ and 
a required SFR of SdOtlf Mq yr^^. 

The SFR values are summarized in Table |2| The SFR re- 
quired by CCSNe is consistent given uncertainties with the 
cosmic SFR extrapolation and the non-rotating UV-SFR, but 
does not overlap within the uncertainties with all other SFR 
estimates. However, the uncertainties are large because of un- 
certain dust corrections on the SFR. If SN 2008S-like events 
are included, there is an even stronger tension including with 
the non-rotating UV-SFR. The Ha-SFR show a more signifi- 
cant tension than the UV, but as we have discussed, the CCSN 
rate probes SFR on time scales that are longer than the Ha 
indicator and more comparable to the UV. 

These SFR values required by CCSNe can be reduced by 
adopting a larger mass range for CCSNe. However, even with 
extreme values of AI„iin — 7Mq and Mmax = IOOM0, the 
required SFR is about 2OOM0yr-i (28OM0yr-i) exclud- 
ing (including) SN 2008S-like events. One possible solution 
is to furth er decrease the minimu m mass for CCSN. This was 
studied bv lBotticella et alj (1201 2b . who suggested the CCSNe 
and non-rotating Ha-SFR are consistent if Mnun = 6Mq. 
Although a similarly small Mnun would be compatible with 
the rotating predictions, we do not favor this scenario because 
-^Aiiin has been accurately determined using progenitor obser- 
vations of CCS Ne in the nearby ^ 20 Mpc volume, i.e., a 
similar volume (ISmartt et al.ll2009b . 

In Figure |5] we show the build up (cumulative) CC- 
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Fig. 5. — The build up of CCSNe (i.e., cumulative), showing the discovered 
CCSNe with (dotted) and without (dot-dot-dashed) SN 2008S-like events, 
predictions from the SFR for Ha (thin solid) and UV (thick solid) for non- 
rotating (black) and rotating (red) calibrations, and prediction based on the B- 
band (dashed), as labeled. All discoveries and predictions have been subject 
to the LVL angular cuts and corrected upwards for the missing volume (see 
text). The total and uncertainty within 11 Mpc are shown by the data point 
with error bars on the right hand side. There are more CCSNe discovered 
than predicted, in particular if the rotating calibration is assumed. 

SNe as a function of distance, and the value within 11 
Mpc as data points. We show the observed CCSNe with 
and without SN 2008S-like events, and predictions based 
on the various SFR estimates, as labeled. We also show 
another CCSN rate estimate based on the galaxy abso- 
lute B-band luminosities using the lat est SNu (1 SN u = 
1 CCSN /IQIOLq^s /century) of LOSS (iLi et al.ll201 ih . The 
SNu depends on the galaxy type as well as the galaxy i?-band 
magnitude. The latter is the rate-size relation found by the 
LOSS group, which introduces a weighting such that smaller 
galaxies have higher SNu (the LOSS rates without the rate- 
siz e relation are in good agreement with previous estimates 
bv lCappellaro et al We use the same galaxy sample as 

we did for the SFR calculation. We correct the galaxy _B-band 
luminosity for foreground Galactic extinction but not for host 
extinction, since the LOSS SNu have not been corrected for 
host extinction either We then sum the CCSN rate from indi- 
vidual galaxies and apply a sky-coverage correction similarly 
to our SFR calculation. 

From Figure |5] it is clear that more CCSNe are discovered 
than predicted. Of the SFR estimates, only the non-rotating 
UV-SFR is marginally compatible with the observed CCSNe 
without SN 2008S-like events. Also, Figure |5] shows that 
the UV-SFR and B-band predictions are quantitatively very 
similar If this is interpreted as the robustness of CCSN rate 
predictions, it leads to the conclusion that there has been a 
chance occurrence of excess CCSNe around ^ 4 Mpc where 
the discovered rate first overshoots the predictions. The statis- 
tics are low in this range; there are only 4 CCSNe within 4 



Mpc. However, comparing this to the i?-band prediction of 
0.85 CCSNe gives a chance occurrence of 0.9% (for the entire 
IIHUGS sample, the respective numbers are 6 CCSNe dis- 
covered and 1.3 expected, a 0.2% occurrence). So, it implies 
a somewhat rare occurrence. Also, it should be noted that the 
cosmic CCSN rate predicted based on the S-band rates is a 
factor ^ 2 smaller than the cosmic SFR rate. Although the 
local and cosmic rates cannot be directly compared, it leaves 
the possibility that the prediction in the local volume based on 
the B-band is also slightly low. 

5. DISCUSSIONS AND SUMMARY 

Based on the new no n-rotational (NOROT) and rotational 
(ROT) stellar tracks of lEkstrom etall (1201 2|) . we used the 
stellar evolutionary synthesis code PEGASE.2 to derive SFR 
calibration factors for non-rotating and rotating stellar popu- 
lations. The ROT population shows significantly increased 
photon output compared to the NOROT population, which 
leads to a decrease in SFR calibration factors of 30% and 
40% for the UV and Ha indicators, respectively (Table [T]). 
C ompared to the w idely-used non-rotating SFR calibrations 
of lKennicutll ( Il998h . the ROT caUbrations are 10% and 40% 
smaller, respectively. These changes to the calibration factors 
are comparable to (or larger than) the uncertainty on the best- 
fit cosmic SFR of HB06 that is 20-30% (l-3cr) at redshifts 
between < z < 1. Thus it is evident that systematic effects 
such as the SFR calibration factors will need to be understood 
better for future precision SFR estimates. 

We compare the SFR estimates to the measured CCSN rates 
to investigate whether the rotational calibrations fit with data. 
At cosmic distances, the best-fit cosmic SFR is too high com- 
pared to the measured CCSN rate data by a factor ^ 2. Stel- 
lar rotation will decrease the cosmic SFR but not enough 
to bridge the difference completely. Furthermore, a part of 
the SFR-CCSN normalization discrepancy may be caused 
by CCSNe that are missed due to a combin ation of intrinsi- 
cally dim and heavily dust o bscured CCSNe dMannucci et al.l 
l2007t iHoriuchi et al.ll20Tll) : this explanation has been sup- 
ported at least at hi gh z by recent CCSN rate measurements 
( iDahlen etaLll2012h . 

At distances of ~ 10 Mpc, several s tudies have found 
hints of an excess of CCSNe (e.g , lAndoetal.1 120051: 



iKistler et alj|201 ll: iHoriuchi et alj|20ru iKistler et al 



l2012l) . 



while iBotticella et al.l ( 120 12h performed an in-depth investi- 
gation and found that the UV-SFR is consistent with the ob- 
served CCSNe. We make new estimates of the SFR and 
CCSN rate in the local 11 Mpc volume, update the dis- 
covered CCSNe, and by comparing the two, find support 
of a n excess of discovered CCSNe (Figure |5]l. Consistent 
with IBotticella et"an ( |2012|) . we find that the total UV-SFR 
is marginally consistent with the observed CCSNe within un- 
certainties. However, upon closer inspection, we find that 
this requires a chance excess of CCSNe in the 4 Mpc vol- 
ume which is rather unlikely (a 0.2-0.9% chance probability). 
The Ha-SFR underpredicts the observed CCSN rate consid- 
erably. Stellar rotation, which decreases both the UV- and 
Ha-SFR estimates, exacerbates the situation and is therefore 
disfavored by current data in the local volume. We empha- 
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size that this result is strengthened by the fact that the locally 
discovered CCSN counts can only increase as surveys probe 
more complete galaxy lists at greater cadence and sensitivi- 
ties. The result is also robust to local overdensities, since both 
the SFR and CCSN rates will be similarly affected. 

Thus, there is a dichotomy between the cosmic CCSN 
rate predicted from the SFR, which overestimates the cos- 
mic CCSN measurements, and the local CCSN rate predicted 
from the observed SFR, which slightly underpredicts the lo- 
cal CCSN observations. Unless the SFR calibrations system- 
atically change with distance, the emerging picture is that 
the SFR calibration factors are constrained to be close to or 
only slightly larger than the currently used values. If it were 
smaller, it would severely underpredict the local CCSN rate. 
If it were larger, it would severely overpredict the cosmic 
CCSN rate. However, it may be slightly larger, because not all 
massive stars need to explode as optically bright CCSNe, al- 
though the fr action of dark collapse s is expected to not exceed 
~ 50% (e.g.. lKochanek et alj|2008l) . 

More generally, the local CCSNe discoveries set a lower 
bound on the SFR, and using this fact we disfavor any signif- 
icant change to the total energy output among massive stars, 
including, but not limited to, rotation. An example is the im- 
pact of binaries. PEGASE.2 treats binarity only for the pur- 
poses of determining chemical enrichment by Type la super- 
novae, but in reality binary interactions affect various aspects 
of stellar evolution, including mixing, mass loss, and surface 
rotation. The binary fraction of massive stars is known to be 
at lea st in the several tens of percent range (e.g., iMason et alj 
and a significant fraction of binaries are estimated to be 
in triple or higher order systems (iGarcfa AMemiiniod'2001'), 
which can lead to interesting inter actions (e.g., . Thompsoni 
I2OIII) . Recentlv. lZhang et al.l(l2012l) found that the UV output 
of stars is increased, leading to a reduction of the SFR calibra- 
tions by 0.2 dex when binaries are included. Such changes 
are also constrained by the observed CCSNe. 

Explanations of the CCSN-SFR dichotomy is itself of in- 
terest. Generally, effects that affect one distance regime and 
not the other are required. Given the simple counting and 
lower limit nature of the locally discovered CCSNe, the most 
obvious option in the local volume is to increase the SFR. 
However, increasing the SFR calibrations will necessarily in- 
crease the cosmic SFR as well. This leaves galaxy dust cor- 
rection as the most realistic explanation. Another possibility 
is that the local galaxy distances are systematically underes- 
timated. Since the SFR scales with the luminosity, and the 
CCSNe are simple counting, a ^ 30% systematic increase 
in the galaxy distances will increase the SFR to match the 
discovered CCSN counts. For the cosmic situation, the cos- 
mic CCSN rate may be too low. One possibility is that the 
cosmic CCSN measurements are mis sing dim CCSNe as ex- 
plained above. iHoriuchi et al.l (1201 ll) shows that the fraction 
of dim CCSNe observed in the local volume is higher than at 
cosmic distances, and argues that incorporating dim CCSNe 
increases the cosmic CCSN measurements to match the cos- 
mic SFR. Alternatively, the cosmic luminosity density used 
to derive many cosmic CCSN rates may be too low, although 
such a systematic change over a wide redshift range is not 



independently motivated. 

Changes to the SFR will cause other effects not dis- 
cussed in this paper For example, the y will affect predic- 
tions of the stellar mas s density (e.g., Madau etal. | 19981 



Madau & PozzettH l2000t ICharv & Elbazl T200ll: ICole et all 
200 ll) . Current predictions based on the HB06 cosmic SFR 
are higher by a factor ~ 2 than the direct measurements 
(iHopk ms & Beacomll2006l) . a fact that has been investigated 
in the context of IMF shape variations and IMF evolution 
dWilkins et al.1 l2008albl) . Rotation decreases the SFR and 
hence the stellar mass density prediction, but the comparison 
to data is more complex. This is because the stellar mass mea- 
surements are derived using galaxy SED fittings, and thus are 
also affected by stellar rotation. Since rotation increases the 
luminosity output of stars, it is likely to reduce the stellar mass 
density measurements also. Therefore, variations of the IMF 
may still be needed. However, quantitative predictions are 
beyond the scope of this paper since it requires calculating 
longer times including effects of low and intermediate-mass 
stars. The SFR h as also been c ompared to the extragalactic 
background Ught dHoriuchi et a l. 2009; Raue & Meyer 2012). 
However, this is not a strong probe of the SFR calibration vari- 
ations, because a smaller calibration is effectively canceled by 
a larger energy output per stellar population. 

The main uncertainty of the present study is the small num- 
ber of CCSNe (between 11 and 22), and the small num- 
ber of star-forming galaxies (the largest 10 galaxies con- 
tribute 45% of the total SFR within 11 Mpc) that make 
dust correction uncertain. Performing comparisons at larger 
distances with more CCSNe and galaxies will dramatically 
improve the situation since even an increase in factor of 2 
in distance would give a factor ^ 8 increase in the rate. 
Presently running observation programs such as the All-Sky 
Aut omated Survey for the Brightest Supernovae (ASAS-SN; 
e.g.. lKhanetani201ll) will find CCSNe in a volume-limited 
sample of nearby galaxies, and the Palomar Transient Fac- 
tory will collect larger numbers of C CSNe in survey mode 
dLaw et alj|2009t IE len &Fieldsll2009l) . Issues such as com- 
pleteness should therefore improve dramatically in the near 
future. Uncertainties inherent to the SFR-CCSN comparison, 
e.g., the proge nitor mass range for CCSN, is rapidly reducing 
(ISmarttll2009l) . CCSNe will provide an excellent and com- 
pelling measure of the SFR in the future that would enable a 
new way of studying various SFR systematic effects. 
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